Vibrational Bloch-Siegert effect in trapped ions 
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When trapped atoms are illuminated by weak lasers, off-resonant transitions cause shifts in the 
frequencies of the vibrational-sideband resonances. These frequency shifts may be understood in 
terms of Stark-shifts of the individual levels or, as proposed here, as a vibrational Bloch-Siegert 
shift, an effect closely related to the usual (radio-frequency or optical) Bloch-Siegert shift and 
associated with rapidly oscillating terms when the Rotating Wave Approximation is not made. 
Explicit analytic expressions are derived and compared to numerical results, and the similarities 
and differences between the usual and the vibrational Bloch-Siegert shifts are also spelled out. 

PACS numbers; 03.75.Be, 32.80.Jz, 37.10.Ty, 37.10.Vz 



I. INTRODUCTION 



A. Bloch-Siegert shift 



The Rotating Wave Approximation (RWA) is applied 
to describe atoms interacting with near-resonant fields. 
It consists on neglecting the rapidly oscillating counter- 
rotating terms in a Hamiltonian. In 1940 Bloch and 
Siegert, studying magnetic resonances, showed that if the 
RWA is not applied, the counter rotating, fast oscillating 
terms give rise to a shift in the resonance frequency of 
the magnetic dipoles, i.e., the Bloch-Siegert shift 
Similar shifts exist in principle in the optical domain, al- 
though much smaller and difficult to detect. 

In the treatment of laser driven trapped ions, apart 
from an optical RWA for the driven optical transition, 
a second or "vibrational RWA" is usually applied [1]. 
Within this approximation, the absorption spectrum of a 
harmonically trapped (two- level) ion consists of a carrier 
band centered at the transition frequency ujq and side- 
bands separated from the carrier by multiples of the trap 
frequency lot- If this second RWA is not applied and vi- 
brational counter-rotating terms are taken into account, 
the energy levels are distorted and thus the position of the 
sideband resonances are shifted, an effect that we shall 
call a "vibrational Bloch-Siegert shift". Its analysis is 
the objective of the present paper in which we shall pro- 
vide explicit expressions within a perturbative approach 
based on the resolvent method. The peculiarities of the 
Hamiltonian for the trapped ion lead to differences be- 
tween the ordinary (in magnetic or optical resonances) 
and vibrational BS effects that we shall also discuss. 

These sideband frequency shifts are important for op- 
timizing the operation of quantum gates based on laser- 



driven trapped ions and can be observed and compen- 
sated experimentally @]- In this context they have 
been understood as the result of Stark-shifts of the in- 
dividual levels due to off-resonant transitions. The con- 
nection between our compact treatment of the frequency 
shift (as a vibrational Bloch-Siegert effect) and the Stark- 
shifts will also be spelled out. Thanks to the system- 
atic treatment with the resolvent method we find gen- 
eral expressions for the frequency shifts of arbitrary side- 
bands and correction terms that had been previously 
overlooked. 



B. Vibrational Rotating Wave Approximation 

We consider a two-level ion with (internal) ground and 
excited states \g) and |e), and transition frequency loq 
among them, moving in an effective one-dimensional har- 
monic potential of frequency lot in x-direction (motion 
in y and z directions is unexcited and ignored through- 
out); this ion is illuminated by a (classical) laser beam 
with traveling wave of wavevector in x-direction and 
wavenumber k^. The system is described in the dipole 
approximation by 

H{t) = Ht + HA + nnB.{<J+ + <J-) cos {uJLt - kLx) , (1) 

with 



rj 



(2) 
(3) 
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where cr^ = |e)(e| - \g){g\, cr+ = \e){g\, and cr_ = |g)(e|; 

is the on-resonance Rabi frequency, which plays the 
role of an atom- field coupling constant, and (a) are 
the creation (annihilation) operators of the harmonic po- 
tential. 

In an interaction picture defined by Ht + Ha and ne- 
glecting the fast oscillating terms within the usual (opti- 



2 



cal) RWA, the Hamiltonian becomes 



HTA{t) 



(*)+i^(t)]g-iAt 



i/.c. , (4) 



where H.c. means "Hermitean conjugate", a{t) 



ae 



and A = ujl—ujq is the detuning 



between the laser frequency and the internal atomic tran- 
sition frequency. The parameter rj = UlXq is known as 
the Lamb-Dicke (LD) parameter, where Xq — yjfi/2mijjj' 
is the extension (square root of the variance) of the ion's 
ground state, i. e., a; = a;o(a + a^). The LD parameter 
is a measure of the trap width on the scale of the laser 
wavelength. Let us denote by \g,n) (|e,n)) the state of 
the ion in the ground (excited) internal state and in the 
n-th motional level of the harmonic oscillator. Using the 
so called BCH identity ff\, 



(5) 



and expanding the exponentials in power series of rj, we 
end up with a Hamiltonian with terms containing a com- 



bination of a±, with 



rotating with a frequency (n — n')ujT 



operators and n' a-operators 



HTAit) 



nn 



2 

H.c. 



_nl ^ {irf) 



E 



n+n 



In/ 1 



n\n' 



(6) 



which in principle couples all the different vibrational lev- 
els via de a^" and a" operators. The combinations sat- 
isfying the A « kuoT {k = n — n') condition will be reso- 
nant, coupling the states \g, n) <-> |e, n+k), while the rest 
of rapidly oscillating terms are usually neglected in a sec- 
ond apphcation of the RWA, a vibrational RWA (VRWA) 
Within this approximation, only co- rotating states 
\g,n) <-> \e,n + k) are coupled and the system becomes 
two-dimensional at each of these resonances. They are 
classified as: (i) A ^ 0: Carrier resonance. The sys- 
tem is described by a Hamiltonian which is equivalent 
to a two level atom in a resonant field, coupling the 
states \g,n) <-> |e,n); (ii) A ^ —kujT- kth red sideband. 
The system is described by a fc-photon Jaynes-Cummings 
type Hamiltonian, which couples |g, n) \e,n ^ k); (iii) 
A ~ kujT- kth blue sideband. The system is described 
by a fc-photon anti Jaynes-Cummings type Hamiltonian, 
coupling the states Ig^n) ^ \e,n + k). 

The differences and similarities of both approximations 
(RWA and VRWA) are discussed in detail in Section 

mm 



II. VIBRATIONAL BLOCH-SIEGERT SHIFT 

If the VRWA is not applied and vibrational counter- 
rotating terms in the Hamiltonian ^ are taken into ac- 
count, a shift in the apparent position of each sideband 




FIG. 1: Bare {Q.r = 0, dashed line) and dressed {Q,r/u)t = 
0.3, solid line) energy levels (in arbitrary units) as a function 
of the laser detuning, with En = nfiLOr being the energy levels 
of the trapping potential. A not too small LD parameter 
r] — 0.4 has been intentionally chosen in order to highlight 
the higher order avoided crossings. 



resonance is observed, an effect that may be compared 
to the standard Bloch-Siegert shift in magnetic or opti- 
cal resonances. This shift is more easily visualized and 
calculated if one writes the original Hamiltonian ([T]) in 
a frame rotating with the field frequency, i.e., in a field- 
adapted interaction picture defined by the zeroth order 
Hamiltonian ^LULaz- In this frame, and after applying 
the usual (optical) RWA, the Hamiltonian becomes time- 
independent. 



H-LA — Tiu'tO' a —a- 



H.C. 



(7) 

The energy levels of this Hamiltonian are plotted in Fig. 
[T]as a function of the detuning. The bare (fifl = 0) levels 
are given by straight lines 



E. 



— ntkoi 



hA 
hA 



(8) 
(9) 



which are degenerate at each of the resonances mentioned 
before (dashed lines in Fig. [1]), but the degeneracies are 
removed and become avoided crossings when the laser is 
turned on (solid lines in Fig. [1]). These anti-crossings 
are nothing but the different sidebands mentioned in the 
previous section and lead to transitions between the dif- 
ferent bare states. In this time independent approach, 
the VRWA corresponds to restricting the theory to a 
subspace spanned only by the two states involved in a 
given anti-crossing neglecting the rest of states, provided 
that the energy splitting of the avoided-crossing is well 
isolated M. 
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FIG. 2: Detail of the energy level diagram at the first blue 
sideband between states \g,0) and |e, I). The bare {Qr = 0) 
energy levels cross each other since the involved states are 
not coupled by the laser (thin-dashed lines) . When the laser is 
turned on (^r/ujt ~ 0.3) and the VRWA applied, the dressed 
energy levels form an avoided crossing, but the position of the 
resonance remains unchanged at Ao = ujt (thin-solid line). If 
the VRWA is not made and counter rotating terms are kept, 
the position of the resonance is shifted to A = Aq -|- Su> (thick 
solid line). In these cases a LD parameter of ?; = 0.1 has been 
used. If the laser is kept turned on but 77 = (the energy levels 
of the harmonic potential do not couple), the states \g,0) and 
|e, 1) are decoupled to all orders in the perturbation. In this 
case, the energy levels cross each other while still being shifted 
by the off-resonant carrier transition (thick-dashed line). 



Nevertheless, a careful treatment of the kth sideband 
(anti-crossings) in which all bare states are included, re- 
veals that the frequency of the resonance does not ex- 
actly coincide with the predicted value A — kujT, but 
is shifted by Sto due to the presence of non-resonant vi- 
brational terms, see Fig. [21 In order to account for this 
effect one has to use a theory which includes not only 
the two main energy levels involved in each resonance, 
but also non-resonant terms. We shall use the resolvent 
method here as described in ■ 



A. Perturbative treatment 

In the low intensity field approximation (flji <^ ^t), 
the Hamiltonian ([7]) describing the system may be writ- 
ten as a sum of a bare (unperturbed) part Hb and a 
(small) perturbative term V , 



Hb 
V{^r) 



,iri{a+a'' ) 



a+ + H.c. 



(10) 
(11) 



The bare states (eigenstates of Hb) of the system satisfy 
HB\a,n)=Ei'^2\c^,n), (12) 

where the index a — g,e accounts for the internal atomic 
state and the bare energy levels Ea n are given in Eq. 
®- 

Let us now consider the \g, Ug) ^ |e, Ue) sideband tran- 
sition, which corresponds to the crossing between the 



bare levels Eg^J^ and Ee%i at the point 
Ao = (ne - ng)ujT 



(13) 
(14) 



of the {E, A) plane, see Fig. [TJ If these levels are close 
to each other but far from other levels, the time evolu- 
tion of the system in the subspace spanned by the states 
\g, Ug) and |e, Ue) may be approximately described by an 
approximate Hamiltonian [9i] 



H = 



Eg^2g + Rgg{Et)) RgeiEo) 

Reg{Ea) Ei%{ + ReeiEa) 



(15) 



where Rap^Eo) = {a,na\R{EQ)\(3,np) are the matrix el- 
ements of the level shift operator defined by 



R{Eo) = PVP 



n=l 



PV 



Q 



Eo-Hb 



V P, (16) 



with P = \g,ng){g,ng\ + |e, 7ie)(e, Jiej and Q = 1 — P. 
Even though it is 2-dimensional, the approximate Hamil- 
tonian (|15|) contains information about all the non- 
resonant states via the operator Q, the projector onto 
the counter-rotating subspace. The eigenvalues of H will 
then give, for weak fields, the correct energy levels, in- 
cluding non-resonant effects.^ It can be easily shown^ 
that the maximum (or minimum) of these perturbed en- 
ergy levels are shifted from their non-perturbed position 
Ao by 



Slu = [R,,iEo) - RggiEo)]/n, 



(17) 



which is the vibrational Bloch-Siegert shift of the position 
of the resonances due to the effect of counter rotating 
(vibrational) terms, see Fig. [51 

Moreover, if the states are coupled by Rge, their corre- 
sponding levels form an avoided crossing with an energy 
splitting given by 



5e=\RgeiEo)\ 



(18) 



^ The resolvent method is capable of providing the exact levels 
substituting Eg by the unknown eigenvalue E in the approximate 
Hamiltonian. Finding E requires iterative procedures, see e.g. 
[lol |. The simplest, explicit, rather than implicit, approximate 
treatment followed here provides the correction to the energy in 
leading order in the perturbation. 

^ Solve dEa^ua, / d./^ = 0, where Ea,na ^'re the eigenvalues of H. 
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where n„,„' = rifl;(n|e"'(°+°^) are the couphng 
strengths between the different vibrational levels of the 
trap [HI [141, see the exphcit expressions in Appendix [Xl 
As we have already pointed out, the approximate Hamil- 
tonian (|15p will only be valid if the anti-crossing is well 
isolated. A criterion for resonance isolation is 6e <^ lut, 
since lut is the energy difference between consecutive res- 
onances. As low intensity lasers {flfj <C ujt) are being as- 
sumed, this criterion will be readily fulfilled and all the 
crossings will be well isolated. 

In order to give an explicit expression of Slo we need the 
matrix elements Rgg and i?ee of the level shift operator. 
To lowest order in the perturbation, they are given by 
(see Appendix \X\ for the explicit calculation) 



Rgg{Eo) 



{g,ng\R\g,ng) = J2 ^ ^(B) ' 
fc=o Eo - h^ ,. 



R\e,ne) = ^ 



iB) ■ (19) 



fc=o Eo- Eg ,, 

k^Tlg 



Then, using the expressions of the bare energies ^ and 
the crossing point of these bare levels (|13p one has from 
P?)) that the position of the \g,ng) <-> |e,rte) resonance 
is given by A = Aq + 6uj, where 



/ 



Slo ~ 



1 

Aujt 



E 

fc=0 



E 



A 



(20) 



which is in principle valid for all values of t] provided low 
intensity fields are assumed. We have verified the va- 
lidity of this expression by comparing with the frequency 
shifts that result from diagonalizing the time independent 
Hamiltonian ([7]) with a large basis of bare states (n ^ 1) 
and numerically finding the maximum or minimum of a 
given energy level, see Fig. [3] (solid lines). 

Note that Eq. (PO)) can also be understood as the sum 
of the Stark shifts for the excited state minus the Stark 
shifts for the ground state due to all non-resonant tran- 
sitions, which provides a connection with previous treat- 
ments [5l,l6|. 



B. LD regime: 77 ^ 1 

A particularly interesting and common regime is the 
so-called LD regime, in which the recoil frequency of 
the ion is much smaller than the trapping frequency, i.e, 
77 ^ 1. Up to quadratic terms in 77, only the same or 
consecutive vibrational levels are coupled, reducing the 
infinite sum in (j20p to one where only "off-resonant car- 
rier terms" (fc = rie in Ree and fc = rig in Rgg) or coupling 
between adjacent vibrational levels is considered. Thus, 




- b) 


1 


1 1 


- n=l 

- /=0 


1 


N _ 
N 


T (1.2 


0.4 


0.6 1) 



n =1 



FIG. 3: Vibrational Bloch-Siegert Shift in units of the trap 
frequency as a function of the LD parameter rj for different 
|(7,ng) |e,ne) transitions and a ratio ^h/uit = 0.01. The 
solid line is a virtually exact result obtained by diagonalizing 
the full Hamiltonian (with a large basis of bare states) and 
numerically finding the maximum or minimum of the corre- 
sponding energy level. This line is indistinguishable from the 
one obtained by computing expression (|20|l to all orders in r). 
The dashed lines correspond to the simplified expression in 
the LD regime, Eq. (|21| . The dotted line in (b) corresponds 
to the shift derived iriQ] and [l^ . 



we find that, see Fig. [3] (dashed lines) 

[l-rf{ng+ 71^ + 1)] 



5uj 



02 



2{ng - ne)ujT 



E 



fc=±i 



Ug + ne + I 
n„ — Ue + k 



0{v'), (21) 



which is valid for Ug ^ rie, since it follows directly from 
(PO)) that the position of the central (carrier, Ug = Ue) res- 
onances are not shifted to any order in rj: the shifts from 
both sides (positive and negative) detunings are com- 
pensated and canceled. (Contrast this with [14| . where a 
different type of shift of the carrier was studied, defined 
by the excitation peak position rather than by the level 
structure). 

The first line in Eq. ipT|) is nothing but the (off- 
resonant) carrier contribution to the vibrational Bloch- 
Siegert shift, while the second line is the contribution 
of the adjacent sideband transitions. It is remarkable 
that in the LD-regime the main contribution to any side- 
band frequency shift is always given by the coupling to 
off-resonant carrier transitions, no matter how far the 
sideband is from the carrier. 

Note that in the strict LD (ry — 0) limit, the states 
\g,ng) and |e,ne) are not connected by the perturbation 
V, and Rge will be zero to all orders. The energy lev- 
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FIG. 4: Bloch-Siegert shift for the first few sidebands, with ex- 
perimental data taken from (JH (see text). The grey scale in- 
dicates the different vibrational levels involved in a particular 
sideband. Black, dark grey, light grey and white correspond, 
respectively, to n = 0, 1, 2, 3. n — Ug for blue sidebands and 
n — He for red sidebands. Note the symmetry between blue 
and red sidebands. The shift is zero for the carrier. 



els will cross while being shifted due to the effect of the 
carrier transitions, which are the only surviving contri- 
butions if f] = 0. The two levels are Stark shifted in 
opposite directions by Rgg — —Ree — ?if^|f/4Ao, shifting 

the position of the resonance by — see Fig. [2] (thick 
dashed Une). 

The vibrational Bloch-Siegcrt shifts can be seen ex- 
perimentally and are indeed quite relevant. In experi- 
ments with single trapped ions to implement quantum 
gates d, H, [l^, ground state sideband cooling [l^ or 
quantum state engineering 0, [13, UBl i the laser has to be 
tuned precisely to a given sideband, so the shift has to be 
taken into account. For recent experiments with ^°Ca+ 
ions [11] {wT = 27r X 1.36 MHz, ^.r = 2t: x 53 kHz and 
rj ~ 0.083), the displacement of the resonant frequency 
at the first blue sideband is as high as one kHz, which 
corresponds to a relative displacement of about 10"^, see 



Fig. H 



nary Bloch-Siegert shift. Compared to previous results 
in [l^ , our formulation is more general and applies to 
arbitrary sidebands. In Q , Steane et al describe the shift 
on the first red sideband as the result of the light shifts 
of the levels because of off-resonant transitions. Aniello 
et al jTsj l , after a series of unitary transformations of the 
Hamiltonian end up with a resonance condition which 
gives a shift also for the first red sideband. In both cases 
the calculated shift is given by 



Slo 



^2 

2ujt 



(22) 



{ug — l,ne = 0), which differs from our result in (PT|) 
by the 77^ correction from off-resonant carrier transitions. 
This correction is indeed easy to miss (in particular, it 
is overlooked if (n|e"'*^°+" ^|ri) is approximated by keep- 
ing only terms up to linear order in rj), but it becomes 
quite significant at moderate values of 77, as shown by the 
numerical comparison of the exact result (from diagonal- 
ization of the Hamiltonian with a large basis), and the 
shifts given by Eqs. ^ and (HH) in Fig. Mh). 



B. RWA vs. vibrational RWA 

In order to highlight the differences between the usual 
RWA and the VRWA (and thus the differences between 
the usual Bloch-Siegert shift and its vibrational version), 
we shall now compare our trapped ion system with the 
interaction of a two-level atom and a quantized field mode 
of frequency co. These two systems are in many ways 
similar if the ion is assumed to be trapped within the 
LD regime [13, [2l|. Since the RWA has to do with co- 
rotating and counter-rotating terms in a time dependent 
Hamiltonian, the different nature of both RWA's is more 
clearly understood if the Hamiltonians describing both 
systems are written in an interaction picture where the 
free evolution of the amplitudes is removed, i. e., 

Hpit) = ^ [(ae-^-* + ate-t)e-oV+ + ff.c.] , 



H' 



TA.LD 



it) 



[(1 -I- iT^ae-"^-^' + iTja^e"^^') e"^*a+ 



+ H.c], 



(23) 



III. DISCUSSION 



where the subscript F refers to the atom-quantized field 
case, and TA,LD to the trapping case, see Eq.®, in 
the Lamb-Dicke regime.'^ The creation and annihilation 
operators have different meaning depending on the sys- 
tem: they increase or decrease the number of photons of 
a given Fock state in the first case or they add or remove 



A. Comparison with previous works 

We have shown that the effect of the vibrational 
counter-rotating terms is to shift the position of the 
sidebands, an effect that may be related to the ordi- 



^ The simplified Hamiltonian HTA,LD{t) would not give the cor- 
rect T]^ correction to the vibrational Bloch-Siegert shift from off- 
resonant carrier transitions, see Eq. I|21| l. but for the qualitative 
comparison of both systems this approximation is enough. 
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a vibrational quantum in the trapped ion case. Despite 
these conceptual differences, both systems look formally 
very similar. 

If near resonant processes are considered {uj coq) in 
Hp, counter- rotating terms like (T„a and cr+a^ can be ne- 
glected since they oscillate with a frequency 2a;, much 
faster than the co- rotating terms, whose frequency of os- 
cillation is S = uj — loq. Counter-rotating terms can never 
be resonant since both the mode frequency iv and the 
transition frequency luq are positive quantities. The main 
effect of keeping these counter-rotating terms is a shift on 
the apparent resonant frequency of the atom luq, i. e., the 
ordinary Bloch-Siegert shift ^] . 

In the trapped ion case, the role of luq is played by the 
frequency difference or detuning A = ujl—ujo, which may 
happen to be positive (blued detuned) or negative (red 
detuned). Terms like cr+a will be resonant if the laser 
is tuned to the first red sideband. This is similar to the 
atom-field coupling case, since after all the interaction 
is described by a Jaynes-Cummings type Hamiltonian. 
However terms like cr+a^ , which would be counter rotat- 
ing in the atom-field coupling case, are resonant in the 
trapped ion system when the laser is tuned to the first 
blue sideband and the system is approximately described 
by an anti Jaynes-Cummings Hamiltonian. 

Another important difference between these systems 
is that for the atom in the quantized-field, transitions 
where only the internal state of the atom changes (with 
no photon absorption or emission) do not happen, while 
in the trapped ion case these carrier-type transitions give 
in fact the main contribution of the LD expansion, i.e., 
the zeroth order in rj contribution in Eq. ([4]). When 
the laser is tuned to a given blue or red sideband, these 
carrier terms become counter-rotating and are respon- 
sible for the leading order contribution in rj of the vi- 
brational Bloch-Siegert shift, Eq. ((2T|) . associated with 
the Stark shift of the energy levels as we have previ- 
ously pointed out. This also explains a factor of two 
discrepancy in the expressions of the dominant terms: 
Sujbs = ^/f/4t^ for the standard Bloch-Siegert shift 
whereas 1(5^^1= il.'j^/2ujT for the vibrational Bloch-Siegert 
shift of the first red or blue sidebands. At variance with 
the vibrational effect the ordinary Bloch-Siegert shift is 
too small to have been observed for optical transitions be- 
cause of the large frequency in the denominator, whereas 
it is relatively easy to observe in the radio-frequency do- 
main 1231. 
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APPENDIX A: CALCULATION OF MATRIX 
ELEMENTS OF THE LEVEL SHIFT OPERATOR 

In this Appendix the matrix elements of the levels 
shift operator are explicitly calculated to leading order 
in the perturbation V, see Eq. pT|) . We will consider 
the crossing between the \g,ng) and |e,ne) bare states, 
with P = \g,ng){g,ng\ + \e,ne){e,ne\ and Q = I - P. 
The first order term (g,ng\V\g,ng) vanishes, since the 
perturbation does not connect states with the same in- 
ternal atomic state. The next order will be given by 



Rgg{Eo) - {g,ng\V 



Eq — Hj 



-y\9,rig) 



(q, nn\V 

2 ^' Eo-Hb 



E 



2 ^ p 



{g,ng\V\e,k) 



E 



Xn„ ,k 



|^^^rag,/c/2| 



E 



(Al) 



where 




Lir:"i(?7^), (A2) 



("> ) being the lesser (greater) of n and n' and 
the generalized Laguerre functions, 

fc=0 ^ ^ 

In an analogous way, Ree is obtained: 



Reem - 2^ - — 



(A4) 
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